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ABSTRACT 

The growth of the supermassive black holes (BHs) that reside at the centres of most 
galaxies is intertwined with the physical processes that drive the formation of the 
galaxies themselves. The evolution of the relations between the mass of the BH, tobh, 
and the properties of its host therefore represent crucial aspects of the galaxy formation 
process. We use a cosmological simulation, as well as an analytical model, to investigate 
how and why the scaling relations for BHs evolve with cosmic time. We find that a 
simulation that reproduces the observed redshift zero relations between rriBH and the 
properties of its host galaxy, as well as the thermodynamic profiles of the intragroup 
medium, also reproduces the observed evolution in the ratio ttibh/'ti* for massive 
galaxies, although the evolution of the mBH/c relation is in apparent conflict with 
observations. The simulation predicts that the relations between tobh and the binding 
energies of both the galaxy and its dark matter halo do not evolve , while the ratio 
"iBH/whaio increases with redshift. The simple, analytic model of iBooth fc Schavg 
(|2010t ). in which the mass of the BH is controlled by the gravitational binding energy 
of its host halo, quantitatively reproduces the latter two results. Finally, we can explain 
the evolution in the relations between tobh and the mass and binding energy of the 
stellar component of its host galaxy for massive galaxies (m* ^ 10^"'^ M0) at low 
redshift {z < 1) if these galaxies grow primarily through dry mergers. 

Key words: Cosmology: Theory ~ Galaxies: Active - Galaxies: Evolution - Galaxies: 
Formation - Hydrodynamics - Galaxies: Quasars: General 



^ 1 INTRODUCTION 



Over the past decade it has become clear that the su- 
permassive black holes (BHs) found at the centres of 
virtually all galaxies with spheroidal components, have 
masses t hat are coupled to the properties of their host 
Ralajcies JMagorrian etal.l 119981: Ferrarese fc Merritt 200d: 



Trernaine et all I2OO2I : iHaring fc Rixl |2004| : iHopkins et all 
20071 ). Additionally, there exists evidence that BH masses 
are coupled to the pro perties of t he da r k matter haloes 

in which they reside (jFerraresel |2002| : iBooth fc Schavd 

2OI0I). Further correlat ions between quasar activity (e.g. 
Boyle fc TerlevichI 119981 ) and the evolution of the cosmic 
star formation rate (e.g. iMadau et al.l 119961 ) provide evi- 
dence that there exists a link between galactic star formation 
and accretion onto a central AGN. 

It has long bee n recognised that the growth of BHs is 
likely self-regulated (|Silk fc Reeslll99a ) and that these tight 
correlations indicate that the growth of BHs is tightly in- 
tertwined with the physical processes that drive galaxy for- 
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mation. However, despite a wide variety of theoretical and 
observational studies, the origin of these relations is still de- 
bated. The study of the evolution of the BH scaling relations 
therefore represents a crucial aspect of the galaxy formation 
process that may provide us with additional clues regard- 
ing the physical processes that give rise to the BH scaling 
relations. 

Addressing these questions observationally is challeng- 
ing. Due to their extremely high luminosities, bright quasars 
provide a promising route to measuring BH masses at high 
redshift through the widths of low-ionization lines that are 
associated with the broad-line region close to the BH and 
using the assumption of virial equilibrium (e.g. IVestergaardI 
I2OO2!). It has, however, been claimed th at this procedure 
syste matically underestimates BH masses (jjarvis fc McLurd 
[2003) • Measuring galaxy masses for these objects is very dif- 
ficult as the BH outs hines the galaxy by a large factor (see 
e.g. the discussion in iMerloni et al.ll201Cr ). Since AGN sur- 
veys are biased towards more massive black holes, selection 
efi^ects also need to be tak en into account (e.g. Shen fc Kclly| 
l2009l : lBennert et al.ll2010l ). which can make it difficult to dis- 
tinguish between evolution in the normalization and in the 
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scatter in the scaling relations (JLauer et al.ll2007l ). In spite of 
these difficulties, measurements of the BH scaling relations 
hav e been made as far ou t as redshift three. 

iMcLure et al.1 (|2006l ) found that the BHs associated 
with radio loud AGN residing in galaxies of a given stellar 
mass are a factor of f our more mas s ive at redshift two than 
in the local Universe. iDecarli et alj (|2010r ) studied the CIV 
line associated with the quasar broad line region in R-band 
selected hosts at both redshifts zero and three and found 
that BHs are typically a factor of seven more massive at 
high redshift for a given galaxy mass. T hese results are con- 
sistent with other observational studies llWalter et al 



Peng et al.ll2006allbl: [Merloni et al.ll2010l : lGreene et al 
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Bennert et al.ll2010l ). Taken together, these papers suggest 



an emerging consensus that at higher redshift BHs in hosts 
of a given mass are systemati cally more m a ssive than in 
the local Universe, although see lJahnke et al.l ((20091) for one 
study that finds no significant evolution. 

The evolution of the relation between BH mass, tubh, 
and stellar velocity dispersion, a,, has been studied util- 
ising the width of the QUI line as a p roxy for stellar 
velocity dispersion (jNelson fc Whittlelll99a ). These studies 
suggest that the ttibh — o"* relation either does not evolve 
(Shields et al. 2003; Gaskell 2009), or doe s so weakly, with 
BHs ~ 0.1 — 0.3 dex more massive at z = 1 (ISalviander et al.l 
l2006l : iGu et al.|[2009l : IWoo et al.|[20oi : iTreu et aLlbOOTp . 

The evolution of the BH scaling relati ons has also been 
studied using numeric al simulations (e.g. iRobertson et al.l 
20061: 1 Johansson et al.ll2009) and semi-analytic models (e.g. 



Malbon et al. 2007: Lamastra et ai]|2010l:lK"isaka fc Koiimal 



2010f ). [Robertson et al.l ( 20061 ) employed simulations of ide- 
alised galaxy mergers, initialised to have properties typical 
of merger progenitors at various redshifts, to construct the 
relation between galaxy stellar mass, m,,, and a, as a func- 
tion of redshift and found that, at a given value of cr*, the 
corresponding jtibh decreases mil dly with increasing red - 
shift. At z = 1 the simulations of iDi Matteo et al.1 (|2008l ') 
have BHs that lie slightly above the z — normalization 
of the 771BH — o" relation. However, these simulations were 
stopped at z — 1 and so cannot inform us about the evo- 
lution of the TTiBH — cr* toward lower redshift. However, for 
z > 1 they predict a weak evolution in the ttibh — cr* relation 
such that at higher redshift galaxies of a g iven velocity dis- 
persio n contain slightly less massive BHs. [Johansson et al.l 
((2OO9I) employed similar numerical techniques to argue that 
it is unlikely that BHs are able to form significantly be- 
fore their host bulges. Semi-analytic models that reproduce 
many redshift zero properties of galaxies also predict that, 
at a fixed a*, BH m asses decrease with increasing redshift 
([Malbon et al.[ [2003). These theoretical models thus pre- 
dict evolutionary trends that go in the opposite direction 
to those inferred fro m observations. Finally, the models of 
[Hopkins et al.[ ([20091 ') predict that, at a fixed stellar veloc- 
ity dispersion, BH masses at higher redshift are either the 
same (for ttibh ~ 10* M©) or slightly more massive (for 
m-BH > 10* Mq) at fixed u* than their redshift zero coun- 
terparts, in agreement with observation. 

On the other hand, the relation between BH mass and 
galaxy bulge mass shows a positive evolution in both semi- 
analytic models (Malbon et_al. 2007; Hopkins et al. 2009) 
and numerical simulations ([Pi Matteo et al.[[2008l ) , the mag- 
nitude of which is comparable to that observed. The larger 



spread in the predictions for the evolution of the msH — o"* 
relation may reflect that it is more difficult to predict veloc- 
ity dispersions, which depend on both mass and size, than 
it is to predict masses. 

In [Booth fc Schavd (|2009l . hereafter BS09) we pre- 
sented self-consistent, hydrodynamical simulations of the co- 
evolution of the BH and galaxy populations that reproduce 
the redshift zero BH scaling relations. These same simula- 
tions also match group temperature, entropy and metallic- 
ity profiles, as well as the stellar masses and age distribu- 
tions of brightest group galaxies (McCarthv et al. 201(J). In 
[Booth fc Schavd ([20101 ) (hereafter BSIO) we used the same 
simulations, as well as an analytic model, to demonstrate 
that rriBH is determined by the mass of the dark matter 
(DM) halo with a secondary dependence on the halo con- 
centration, of the form that would be expected if the halo 
binding energy were the fundamental property that controls 
the mass of the BH. In the present work we use the same 
models to investigate why and how the BH scaling relations 
evolve for massive galaxies. 

This paper is organised as follows. In Sec. [2] we sum- 
marise the numerical methods employed in this study and 
the simulation analysed. In Sec. [3] we present predictions for 
the evolution of the BH scaling relations and compare them 
to observations. We find that the evolution in the ttibh — rn, 
relation predicted by the simulations is in excellent agree- 
ment with the observations, while the measured weak evo- 
lution in the ttibh — o" relation is in apparent disagreement, 
and predict that while BH mass increases with redshift for 
fixed halo mass, the relations between tttbh and the binding 
energies of both the host galaxies and DM haloes do not 
evolve. We demonstrate in 3] that the a analytic description 
in which ttibh is coupled to the DM halo binding energy 
can reproduce the evolution of the relation between BH and 
halo mass. Furthermore, we show that the evolution in the 
relations between the BH and the stellar mass and binding 
energy can be understood in terms of the more fundamental 
relation with the binding energy of the dark halo and the 
growth of massive galaxies through dry mergers. Finally, we 
summarise our main conclusions in Sec. [5] 



2 NUMERICAL METHOD 

We have carried out a cosmological simulation using a signif- 
icantly extended version of the parallel PMTree-Smoothed 
Particle Hy drodynamics (S PH) code GADGET ill (last de- 
scribed in [Springell [2005h . The simulation and code are 
described in detail in BS09, we provide only a brief sum- 
mary here. In a ddition to hydrodynamic forc es, we treat 
star form ation ([Schave fc Dalla V ccchia 2008!), supernova 
feed back ([Dalla Vecchia fc Schave 2008 ). radiative cool- 
ing |\yiersma et aL[[2009ah . chemodynamics ([Wiersma et al.[ 
2009bl ) and black hole accretion and feedback (BS09, 



Springel et al.| [200a). We summarise in Sec. 12.11 the essen- 
tial features of the BH model. 

The properties of central galaxies and DM haloes 
are calculated by first identifying the most gravitationally 
bound particle in each DM halo using the algorithm SUBFIND 
(Springel et al. 2001; Dolag et al. 2009), which is then con- 
sidered the halo centre. All stars within a radius of O.lSrhaio 
are then assigned to the central galaxy. We note that our 
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conclusions are insensitive to the exact choice for this ra- 
dius, and whether we use a fixed physical value or a fixed 
fraction of the halo virial radius. As long as the sphere en- 
closes the central object, our results are insensitive to this 
choice. Halo mass, rrihaio, is calculated as the total mass 
enclosed within a sphere, centred on the most bound par- 
ticle in the halo, that has a mean density of 200 times the 
mean density of the Universe and the virial radius rhaio is 
the radius of this sphere. Because it is not expected that the 
same physics holds for both the central galaxy of a halo and 
its satellites, which are expected to rapidly have their gas 
supply stripped when they become satellites, our analysis is 
restricted to BHs identified as residing in the central galaxy 
in a DM halo, defined as the galaxy closest to the centre of 
the DM potential well of each halo. 

2.1 The black hole model 

Seed BHs of mass rrisccd ~ 10~^7ng ~ 10^ M0 - where rrig is 
the simulation gas particle mass - are placed into every DM 
halo that contains more than 100 DM particles (which cor- 
responds to a DM halo mass of 4.1 x lO'^" Mq/H) and does 
not already contain a BH particle. Haloes are identified by 
regularly running a friends-of-friends group finder during the 
simulation. After forming, BHs grow by two processes: accre- 
tion of ambient gas and mergers. Gas accretion occurs at the 
minimum of the Eddington rate, mEdd = iTrGmBHrnp/eiaTC 
and rhaccr ~ a'i'rTG'^m^-^p/{cl + v^)^'^ , where rUp is the pro- 
ton mass, ctt is the Thomson cross-section, e^ is the radia- 
tive efficiency of the BH, c is the speed of light, Cs and p 
are the sound speed and gas density of the local medium, v 
is the velocity of the BH relative to the ambient medium, 
and a is a dimensionless efficiency parameter. The parame- 
ter a accounts for the fact that our simulations possess nei- 
ther the necessary resolution nor the physics to accurately 
model accretion onto a BH on small scales. Note that for 
Q = 1 this accretion rate r educes to the so called Bondi- 
Hoyle (|Bondi fc Hovlelll944l l rate. 

As long as we resolve the scales and physics relevant 
to Bondi-Hoyle accretion, we should set a = 1. If a simu- 
lation resolves the Jeans scales in the accreting gas, then it 
will also resolve the scales relevant for Bondi-Hoyle accre- 
tion onto any BH larger than the simulation mass resolution 
(BS09). We therefore generally set a equal to unity. How- 
ever, this argument breaks down in the presence of a multi- 
phase interstellar medium, because our simulations do not 
resolve the properties of the cold, molecular phase, and as 
such the accretion rate may be orders of magnitude higher 
than the Bondi-Hoyle rate predicted by our simulations for 
star-forming gas. We therefore use a power-law scaling of 
the accretion efficiency such that a = (wh/^Ih)'^ in star- 
forming gas, where n^ = 0.1 cm~'^ is the critical density 
for the formation of a cold, star-forming gas phase. The pa- 
rameter /3 is a free parameter in our simulations. We set 
/3 = 2, but note that the results shown here are insensitive 
to changes in this parameter when P > 2. We note that we 
do not resolve the Bondi radius of BHs less massive than 
the particle mass in our simulations, and that for BHs with 
rriBH ~ rng the Bondi radius is unresolved unless the density 
is low or the temperature high. Our choice of a = 1 there- 
fore provides an underestimate of the true accretion rate in 
these regimes. However, even setting a = 100 for all densi- 



ties gives very similar results (jBooth fc Schavell2009l ). This 
is because all BHs accrete almost all of their mass in short, 
(near) Eddington-limited bursts of accretion and thus that 
our treatment of accretion in low-density environments is 
less important. The accretion model in high-density environ- 
ments is necessarily very crude, but we note that our results 
are insensitive to the details of the accretion model as long 
as Q is sufficiently large that the BHs become more massive 
than observed in the absence of feedback and if two reason- 
able conditions are met that are necessary for self-regulation 
to be possible. Firstly, the BH accretion rate must increase 
with increasing density, and secondly it must increase with 
BH mass (see BS09). 

Energy feedback is implemented by allowing BHs to in- 
ject a fixed fraction of the rest mass energy of the gas they 
accrete into the surrounding medium. The energy deposition 
rate is given by 



E 



eftrmaccrC 



EfEr . 2 

; mBHC , 



(1) 



where rhaccr is the rate at which the BH is accreting gas, 
and rriBH is the rate of BH mass growth. 

We set Er to be 0.1, the mean value for ra- 
diatively efficient accret ion onto a Schwarzschild BH 
(jShakura fc Sunvaevlll97 3) and use ef — 0.15 as our fiducial 
value. It was shown in BS09 that, for ef = 0.15, this simu- 
lation reproduces the observed redshift zero rriBH — m-* and 
m-BH — f* relations. Energy is returned to the surroundings of 
the BH 'thermally', by increasing the temperature of A^hcat 
of the BH's neighbouring SPH particles by at least ATmin. A 
BH performs no heating until it has built up enough of an en- 
ergy reservoir to heat by this amount. Imposing a minimum 
temperature increase ensures that the radiative cooling time 
is sufficiently long for the feedback to be effective. We set 
A^'hcat = 1 and ATmin ~ 10* K but the results are insensitive 
to the exact values of these parameters (see BS09). 



2.2 The cosmological simulation 

The simulation employed in the current work uses a cu- 
bic box of size 50 comoving Mpc/h and assumes periodic 
boundary conditions. The simulation contains 256"^ particles 
of both gas and collisionless cold DM and is evolved down to 
redshift zero. The DM and initial baryonic particle masses 
are 4.1 x IO^Mq/Zi and 8.7 x lO^'Mo/fe, respectively. Co- 
moving gravitational softenings are set to 1/25 of the mean 
interparticle separation down to 2 = 2.91, below which they 
switch to a fixed proper scale of 2kpc/ft. The simulation 
employed in this study was previously also analysed as the 
fiducial simulation in BS09 and BSIO 

Comparison of the simulation employed in this study 
to an otherwise identical one with eight (two) times lower 
mass (spatial) resolution informs us in what mass and red- 
shift range the relations between ttibh and galaxy and halo 
properties are numerically converged. The relation between 
TTiBH and both m, and rrihaio is numerically converged up 
to redshift two for all haloes with rriBH > lOmsocd. Mea- 
surements of stellar velocity dispersion are, however, only 
converged for z < 1 and rriBH > lO'^ rUsacd- We will only give 
results for redshifts and BH masses for which the results are 
converged with respect to numerical resolution. 

We note that these simulations do not resolve the scales 
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Figure 1. The median m^n/mt ratio as a function of redshift 
for galaxies of different stellar masses (left panel) and for BHs 
above a given mass (right panel) predicted by the cosmological 
si mulation. The points with error bars show the measurements 
of lMerloni et al.l 1 120101 . red triangles) for a sample of galaxi es in 
the m ass range 10.5 < log(ni,/M0) < 11.5 and Dccarli c t al.l 
1120101 . blue diamonds) for a sample of galaxies with typical stel- 
lar masses ~ 10^^ M©. We show a number of different mass cuts 
for the simulation results to demonstrate that the results are in- 
sensitive to the particular mass cut chosen. Both the simulations 
and observations in this plot show the total stellar mass of the 
galaxies. The simulation predicts that the median m-Qn/mt ratio 
increases with redshift. The predictions are in excellent agreement 
with the observations, regardless of the particular mass cut made. 



on which the BH is the gravitationally dominant compo- 
nent in the galaxy and so cannot be used to study BH self- 
regulation on the smallest scales. However, the simulations 
do have sufhcient resolution for baryons to be gravitationally 
dominant in the centres of haloes. We cannot conclusively 
rule out that if we increased our mass resolution signifi- 
cantly and used more sophisticated sub-grid models that the 
BH woul d self-regulate on different scales. However, sugges- 
tively, in lBooth fc Schavd ([2010) we verified that in simula- 
tions with a spatial resolution of 0.5 kpc/Zi that, at z = 2, the 
BH masses scale in the same way as in the lower-resolution 
simulations. 



3 SIMULATION RESULTS 

3.1 The evolving relations between black holes 
and galaxies 

Fig. [l] compares the predicted evolution of the median 
tibh/jti, ratio for different minimum stellar (left panel) and 
BH (right panel) masses with observations of AGN in galax- 
ies with 771* ~ 10^^ M0. The ttibh/tti* ratio increases with 
redshift, in close agreeme nt with observations l| Merloni et al.l 
l2010l : lDecarli et al.ll2010l ). At redshift zero this agreement is 
unsurprising because the efficiency of AGN feedback in the 
simulation was tuned to reproduce the normalisation of the 
z = Q BH scaling relations. The agreement at higher redshift 
represents, however, a non-trivial prediction of a model in 
which BHs self-regulate their accretion through the coupling 
of a small fraction of the radiative energy to the ambient 
medium. 

In Table [T] we show the predicted evolution up to z = 



Table 1. Evolution of the normalisation of the power-law rela- 
tions between BH mass and the properties of its host. The sample 
includes all galaxies that contain a BH with 7r7BH > 10^ TTis^cd ~ 
10^ M0, corresponding to 162 (132) haloes at 2 = (z = 1). 
The median BH, stellar, and halo mass are 10^'', lO^'^'®, and 
10^^'® Mq, respectively. The central column gives the median 
change in logjg ttibh a-t redshift one relative to 2 = for fixed 
values of the quantity listed in the left column. The right column 
shows the slope, a^, of the best fit power-law describing the rate 
of evolution of the scaling relation (Eq. [2]| over the redshift range 
0—1. From top to bottom, we consider evolution of ttibh for fixed 
stellar mass, halo mass, central stellar velocity dispersion, galaxy 
binding energy (mtcrl), and halo binding energy (Eq. [Sj. Errors 
are calculated from 10^ bootstrap resamplings of the data. 



Variable 


A 


^°glO mBH(^ = 0) 


Qs 


771, 




0.20 ±0.05 


0.52 ±0.05 


nihalo 




0.23 ±0.03 


0.65 ±0.06 


C* 




-0.09 ± 0.04 


-0.32 ±0.05 


U* 




0.02 ± 0.03 


0.05 ±0.06 


f^halo 




0.01 ±0.03 


0.03 ±0.05 



1 of the amplitude of the relations between ttibh and the 
masses, velocity dispersions and binding energies of both the 
host galaxies and the host DM haloes. The central column 
gives logiQr7iBH(z = 1) — logj^Q 777bh(2 = 0), calculated by 
fitting power-law relations under the assumption that the 
slopes of the scaling relations do not evolve, which is a good 
approximation in the redshift range studied herqj. The right- 
most column of Table [l] gives the slope, as, of the power- law 
evolution in the amplitude of each scaling relation 






0c(l + 2)'^ 



(2) 



where X is one of the variables listed in the left column, and 
770 is the slope of the ttibh — X relation at 2; = 0. We find 

relation, in good 



m. 



that Qfs = 0.52 ± 0.05 for the ttib h 

agreement with iDi Matteo et al.l (|2008l ). who found a^ — 
0.5. 

The evolution of the ttibh — o", relation is smaller but 
significant, with as = —0.32 ± 0.05. This is in apparent 
disagreement wit h various observational s tudies that either 
infer a positive (ISalviander et all I2OO6I : iTreu et al.) l2007l: 
Woo et aLlbOoil: IGu et al.ll2009l ) or negUgible (IShields et all 
20031 : iGaskelll [2009|) evolution in the normalisation of the 
tttbh — cr, relation. The predicted e volution does, however , 
agree with other simulatio n studies (Robertson et al.l 12006! ) 
and semi-analytic models (jMalbon et al.l 12007 ). Taken to- 
gether, the simulation predicts that the tttbh — ttt, and 



^ At 2 = (2 = 1) the slope of the relation between ttibh ^^nd 
C/halo is 1.01 ± 0.14 (0.96 ± 0.17) and the slope of the relation be- 
tween TTiBH and U* is 0.93± 0.07 (0.96±0.09), con siste nt with the 
2 = observational results of lFeoli k. Meld ||2007| ') and lFeoli et al.l 
1120101 ) . At all redshifts these slopes are consistent with unity. The 
slopes of the relations between 777bh and 777, , cr* and Trxjiaio are 
1.16±0.06 (1.2±0.2), 4.6±0.8 (4.4±0.8) and 1.5±0.2 (1.5±0.3). 
There is thus no evidenc e for evolution in any of the slopes, which 
agrees with the re sults of|Robertson_et_alJ 1(200(1) for the ttibh ~o"* 
relation and with lDi Matteo et al.l 1 120081 ) for the 777bh — 777, rela- 
tion. T he slopes of the 2 = relations are consistent with obser- 
vatio ns llHaring fc Rixll2004l : iTremaine et al.ll2002l : lBandara et al.l 
l2009t) . 
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rriBH — <y* relations evolve such that the relation between 
JTiBH and stellar binding energy (ex 771*(t^) is independent of 
redshift (a^ = 0.05 ±0.06). 

It is tempting to conclude from the finding that the 
ratio rriBu / i'mt.a^) does not evolve that the BH mass is de- 
termined by the binding energy of the galaxy. However, we 
demonstrated explicitly in BSIO that the BH mass is instead 
controlled by the binding energy of the DM halo. This im- 
plies that the binding energy of the galaxy tracks the binding 
energy of the halo, which we will confirm and explain below. 



3.2 The evolving relations between black holes 
and dark matter haloes 



Model with c(m ,z) 
Model with c(m ,z = 0) 
Simulation 




We now turn our attention to the relations between the mass 
of the BH and the DM halo in which it resides. In BSIO we 
argued that a BH grows until it has injected an amount of 
energy into its surroundings that scales with the binding 
energy of its host DM halo. We therefore do not expect the 
JTiBH — t/haio relation to evolve. Indeed, in the simulation the 
amplitude of this relation is independent of redshift, with 
Os = 0.03 ± 0.05, as would be expected for a fundamental 
link between rriBH and the binding energy of the host DM 
halo. 

We do, however, expect the relation between tubh and 
n^halo to evolve. At higher redshift, haloes of a given mass 
are more compact than their redshift zero counterparts and 
are thus more strongly gravitationally bound. This means 
that, at a fixed halo mass, more energy is required to eject 
gas from haloes at high redshift and in order to self-regulate, 
BHs must grow to be more massive. Fig.[2]shows the normal- 
isation of the rriBH — Jrihaio relation as a function of redshift 
and confirms that the simulation predicts the amplitude of 
this relation to increase with redshift (red, dashed curve), 
with Qs = 0.65 ± 0.06. 



1-Fz 



Figure 2. Predicted evolution of the normalisation of the rriBH — 
mijaio relation. The red, dotted curve shows the median evolution 
predicted by the cosmological simulation when all BHs for which 
both BH and halo properties are well resolved (mBH > 10^ rrigeed) 
are included. The grey region represents the allowed range in evo- 
lution predicted by the analytic model of BSIO, which assumes 
that the BH mass is controlled by the binding energy of the DM 
halo. The binding energy of an NFW halo depends on mass, red- 
shift, concentration (which itself depends on both mass and red- 
shift) and on the radius, r^j, at which it is evaluated. The grey 
region corresponds to rgj/rhaio = 0.1 — 1.0 (bottom to top) and 
the solid black line to r^j = 0.22rijaio) the value for which we 
which we predict rriBH oc "m^^^i^ at 2 = 0, in accord with both 
observations (Bandara et al. 2009) and simulation (BSIO). For 
comparison, the black dashed line shows the evolution that our 
analytic model would have predicted if we had ignored the evo- 
lution of the c(mijaio) relation. At all redshifts the normalisation 
of the simulated mBH ~ 'T*halo relation (red, dotted curve) agrees 
with that predicted by the analytic model based on the assump- 
tion that the fundamental relation is between BH mass and the 
binding energy of the DM halo. 



4 EXPLAINING THE EVOLUTION 

As we already noted, the idea that the binding energy of 
the dark halo controls the mass of the BH explains our find- 
ing that the mBH — f^haio relation is independent of red- 
shift. We will now show that the analytic model of BSIO 
also reproduces the evolution of the rriBH — m-haio relation 
and that it can explain the observed evolution in the scaling 
relations with the stellar properties if the observed galaxies 
evolve predominantly through dry mergers, as predicted by 
the simulation. 



4.1 The mBH 



T^halo 



relation 



If the energy injected by a BH is proportional to the halo 
gravitatio nal binding energy, [ /halo, then, for a DM halo with 
an NFW (|Navarro et al.lll997h density profile (BSIO) 



mBH OC (/halo tX 



Thslo 



OC f{c,x){l + z)m^^'^^^, (3) 



5/3 



where mhaio is the halo mass, c(m,haio, z) is the halo concen- 
tration, X is defined to be a; = rej/rhaio, ''cj is the physical 
scale on which BH self-regulation takes place, and /(c, x) is 



the function 



(ln(l + c)-c/(l + c))' 

1 21n(H-CT) 



(l + cx)2 



1 -I- ca; 



(4) 



Simulations have shown that c is a function of both redshift 
and halo mass, and scales approximately a^J (JDuffv et al.l 
l2008h 



^halo^(l + ^)" 



(5) 



Combining Eqs. |3]-E1 BSIO found the slope of the thbh — 
Tihaio relation to be weakly dependent on Toj. At z = it 
varies from no = 1.50 for x = 0.1 to no = 1.61 for x = 1.0. In 
order to exactly matc h the slope of 1.55 that is both observed 
(|Bandara et al.ll2009l find 1.55 ± 0.31) and predicted by the 



^ Note that the slope of the power-law dependence o f concentra- 
tion on redshift depends on the halo definition used I IDuffv et al.l 
l2008f) . 
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simulations (BSIO find 1.55 ± 0.05 for the same simulation 
as is analyzed herqj), we would need to use x = 0.22. 

By using an NFW density profile and Eq. ((5|, we have 
implicitly assumed that the dark matter profile is well de- 
scribed by the results obt ained from simul ations that in- 
clude only dark matter. iDuffv et al.l l|2010l ) have recently 
shown that, on the scales of interest here, the back-reaction 
of the baryons onto the dark matter is in fact very small 
if feedback from AGN is included, as required to reproduce 
the observed s tellar an d gas properties of sroups of galajcies 
jMcCarthv et al. 201Q: iPuchwein et all l2008l : IPabian et all 
l2010l : lDufrvet al.ll2010l '). 

If, as argued in BSIO, the BH mass is controlled by the 
DM halo binding energy, then we expect the rriBH — »Tihaio 
relation to evolve because the halo binding energy depends 
not only on halo mass, but also on the virial radius and 
concentration, both of which vary with redshift for a fixed 
halo mass. 

If the c — rrihaio relation did not evolve, then we would 
expect rriBH (jTihaio) oc (l + z) (Eq.[3]). However, because halo 
concentration decreases with redshift (Eq. O we expect the 
actual evolution of the niBH — wihaio relation to be weaker, 
i.e. Oa < 1. The resulting relation between BH mass and DM 
halo binding energy predicts that, at a given nihaio, JTIbh in- 
creases with redshift and that by 2: = 2 BHs are between 
1.5 (for rej/rhaio = 0.1) and 2.6 (for rej/rhaio = 1-0) times 
more massive than at redshift zero. For our fiducial radius 
of self-regulation of 2; = 0.22, BHs are 2.1 times more mas- 
sive, in excellent agreement with the simulation prediction 
of as = 0.65 ± 0.06 (Table [T|). 

The evolution predicted by Eqs.[3]-[5]is shown in Fig. [2] 
The grey shaded region outlines the analytic prediction for 
the evolution in BH mass over the range rcj/rhaio ~ 0.1 — 1.0 
and the solid black line shows the prediction for rcj/rhaio = 
0.22 (the value that reproduces the slope of the redshift zero 
m-BH — Tihaio relation). The red, dotted curve shows the sim- 
ulation prediction for the evolution of the rriBH — n^halo re- 
lation, including all BHs with niBH > 10^ m,seed- At all red- 
shifts the normalisation of the simulated ttibh — wihaio rela- 
tion is compatible with that predicted by the analytic model. 
For comparison, the dashed, black line shows the predicted 
evolution of the rraBH — JTihaio relation if c(mhaio) were inde- 
pendent of redshift. The analytic model can only reproduce 
the simulation result if the evolution of the concentration- 
mass relation is taken into account. 

The evolution of the niBH — nihaio relation thus provides 
additional evidence for the idea that the masses of BHs are 
determined by the binding energies of the haloes in which 
they reside. 



4.2 The relations between rriBH and galaxy stellar 
properties 

Considering now only the stellar masses for which the evo- 
lution has been measured observationally (m* ~ 10^^ M0) 
and the redshift range for which all of the stellar and BH 
properties of the galaxies are converged numerically {z < 1), 
we ask if we can explain how the relations between BH mass 



^ We quoted a slope of 1.5 it 0.2. Our error bar is greater because 
BSIO used rriBH > lOnisccd whereas we require rrtBH > lOOrTiseod- 




0.0 0.2 0.4 0.6 O.i 

Redshift 



1.0 



Figure 3. Galaxy properties as a function of redshift for m« ^ 
10^^ M0 (at each redshift we selected the 20 galaxies with stellar 
masses closest to this value). The yellow shaded region shows the 
area that contains the 25th to 75th percentiles of the data and 
the black, solid curve shows the median. The top panel shows the 
ratio of stellar to halo binding energy and the bottom panel shows 
the specific star formation rate (SSFR= mt/m,). The dashed line 
in the bottom panel shows the inverse of the age of the Universe. 
Galaxies with SSFRs below this line may be considered passive. 



and galaxy stellar properties evolve. BSIO showed that the 
BH mass is determined by the binding energy of the DM 
halo, which explains why the mBH — t/haio relation does not 
evolve. We find that the thbh — U, relation also does not 
evolve, implying that, over the range of redshifts and masses 
investigated here, f/, oc C/haio- The top panel of Fig. [3] con- 
firms that this is indeed the case in our simulation. 

For the binding energy of the galaxy to track that of the 
halo, we require the two to grow through the same mech- 
anism. This condition is met if the galaxies grow primarily 
through dry mergers. In the absence of significant in-situ star 
formation, both the stellar component, which is predomi- 
nantly spheroidal for massive galaxies, and the DM halo are 
coUisionless systems and are therefore expected to evolve in 
a similar manner. 

The bottom panel of Fig. [3] shows that at 2; <C 
1 the specific star formation rates (SSFR= rhs/mt) of 
galaxies with m, ~ 10^^ Mq are significantly lower than 
the inverse of the Hubble time, l/iu, implying that 
the galaxies are indeed not growing significantly via in- 
situ star-f ormation, in agreement with various observa- 
t ions ( e.g. 'Schawinski ct al' '2006! : Ivan der Wei et all l2009l : 



201G| ). 



Bezans on et al. 2009; van Dokku m et al.ll201Gl ). The analy- 
sis in Fig. 13] was carried out for galaxies at a fixed stellar 
mass, but the same results hold if we trace individual galax- 
ies through time. While the stellar masses of the most mas- 
sive z — galaxies grew on average by a factor 2.65 since 
2 = 1, the fraction of stars in the redshift zero objects with 
birt h redshifts below z = 1 is, on average, only 15%. 

iMerloni et"al] l|2010t ) studied the evolution of the BH 
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scaling relation in the redshift range 1 < z < 2.2, and found 
that at z = 1 a large fraction of galaxies would be identi- 
fied as star-forming. This is consistent with our results as 
the SSFRs of our galaxy sample are increasing with increas- 
ing redshift so that a significant fraction of them would be 
identified as star-forming at z ~ 1. The fraction would be 
even higher if we had required the galaxies to contain active 
AGN, as is the case for the objects selected in the observa- 
tions, because AGN activity is typically accompanied by a 
temporary increase in the star formation rate. 

One might naively think that the ttibh — wi* relation 
should not evolve at all if the galaxies grew predominantly 
through dry mergers. However, the progenitors of galaxies 
with m» ~ 10^^ Mq at z = typically formed their stars and 
grew their BHs later than the progenitors of galaxies that 
already have the same stellar mass at z = 1. Because the pro- 
genitors of higher redshifts galaxies formed earlier, they have 
higher binding energies and thus greater BH masses relative 
to their halo masses. Thus, even if m« ~ lO'^^ M© galaxies 
are growing predominantly by dry mergers at both redshifts, 
they may have different BH masses. Observe, however, that 
the evolution in the rriBH — rn,. relation that we predict for 
these massive galaxies is only mild (about 0.2 dex; see Fig.[T| 
and that in situ star formation thus become important for 
z> 1 (see Fig. 13). 

If the ratio of the stellar to halo binding energies re- 
mains constant, we can write 



rriBH oc l7haio ccUt oc m.o-, . 



(6) 



To explain the evolution of ttibh /jti, that is observed for 
m* ~ IO^'^Mq, and which the simulation reproduces, we 
need to know how the m* — a, relation evolves for such 
galaxies. 

Measurements of the evolution of the m,* — a, relation 
have so far only been underta ken for a small number of 
objects. ICappellari et al.l (|2009l ) presented stacked observa- 
tions of seven galaxies with m* ~ lO^'^ Mq in the redshift 
range 1.6 < z < 2.0 and found that these galaxies typically 
have the same (t* as the very highest velocity dispersion 
early- type galaxies of the same mass in the local Universe. 
This is in agreement with observations showing that galax- 
ies of a given stellar mass are more compact at higher red- 
shifts (see e.g. IWilliams et al.ll2O10l 'l. In deed, for early-type 
:alaxi es with stellar masses ~ lO^'^ M© ICenarro fc Truiilld 
20091) find that typical velocity dispersions decrease from 
« 240 km/s at 2 = 1.6 to fs 180 km/s at 2 = 0, which 
implies cr* oc (1 -f z)"'^. Combined with Eq. © this yields 
wibh/ot, (X {1 + z)°'^ , in good agreement with the as ~ 0.5 
that is observed for galaxies of this mass. 
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m-BYi/m, oc (1-1-2)"% with a^ ~ 0.5, and m-B-a/crt ot (l-|-2)"° 
with Qs Ri —0.3, in apparent confiict with recent observa- 
tions. The ratio between the BH mass and the binding en- 
ergy of the dark halo is independent of redshift, in agreement 
with BSIO who argued that the BH mass is controlled by the 
halo binding energy. The simple analytic model of BSIO, in 
which the BH mass is assumed to scale in proportion to the 
binding energy of the dark halo, not only reproduces the 
simulated redshift zero msH — fTihaio relation, but also its 
evolution. For a fixed halo mass BHs are more massive at 
higher redshift because the haloes are more compact and 
thus more tightly bound. Assuming an NFW halo density 
profile and the evolution of the halo concentration-mass rela- 
tion predicted by simulations, the model can quantitatively 
account for the predicted evolution. 

The simulation predicts that the ratio between the 
BH mass and the binding energy of the stellar component 
of the galaxy is also independent of redshift (at least for 
m* ^ 10^^ Mq and z < 1), even though BSIO demonstrated 
explicitly that the correlations between BH mass and stel- 
lar properties are not fundamental. This result is, however, 
consistent with a picture in which massive galaxies grow 
primarily through dry mergers at low redshift, which we 
showed to be the case in the simulation. Combined with the 
observed evolution in the m« — ct* relation, this idea can 
quantitatively account for the evolution in the rriBH — m* 
relation. 

One interesting implication of this scenario is that the 
evolution of the relations between BHs and the properties of 
their host galaxies may differ for galaxies that do not grow 
predominantly through dry mergers, as would be expected 
for lower masses and at higher redshifts. We will investigate 
this further in the a future work, employing higher resolution 
simulations. 
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5 CONCLUSIONS 

We have used a self-consistent cosmological simulation that 
reproduces the observed redshift zero relations between m-BH 
and both galaxy and halo properties (as well as the thermo- 
dynamic profiles of the intragroup medium) to investigate 
how, and why, these relations evolve through time. 

The relation between BH mass and host galaxy mass 
predicted by the simulation is consistent with available ob- 
servations at z < 2, which are currently confined to m* ^ 
10^^ Mq galaxies. For such galaxies we predict that the ratio 
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